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Abstract 

We show that the electronic band theory in the momentum space requires information about 
the transport of the eigenfunctions . The transport of the eigenfunctions in the Brillouine Zone 
induces spin connections (gauge field in the momentum space) and curvatures (the equivalent of 
the electromagnetic strength). When the theory is applied to topological materials characterized 
by discrete symmetries such as time reversal,parity inversion, and charge conjugation the curvature 
and the spin connections needs to satisfy constraints conditions. As a result of the constraints the 
curvature generates topological invariants such as the first and second Chern number. The Chern 
numbers are reveald by measuring the response to an external electromagnetic field. 

We will study transport in Topological Insulators,Weyl Semimetals (on flat and curved surfaces) 
and Topological Superconductors for which we compute the Andreev crossed reflection. Using 
a p — wave wire coupled to two rings we show that the persistent current in the rings contains 
information about the Majorana fermions. 
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I-Introduction 


One of the important ideas in Condensed Matter Physics is the concept of topological 
order [Ml EHni HSl [231 IMl ESI EE im EM31 IS2|. Insulators with a single Dirac cone 
which lies in a gap such as Bi 2 Se 3 , Bi 2 Te 3 and Bii^^Sbx represent the experimental 
realization of Topological Insulators {T.I) and recently the Weyl semimetal |1H7| such as 
WTe 2 and HgCr 2 Se 4 . 

At the surface of the three dimensional Topological Insulator {T.I.) one obtains a two 
dimensional metallic surface characterized by an odd number of chiral excitations, due to 
Kramers theorem, electrons are protected against backscattering dl and localization m 
l66] . When time reversal is broken localization effects are observed [65] . The surface physics 
has been realized in CdTe/HgTe/CdTe quantum wells. The quantized spin-Hall effect has 
been proposed [3] and observed m |6T] and recently the Anomalous Hall effect has been 
measured [26]. The spin resolved photoemission [2] has been used to identify the surface 
states. Topological superconductors and their identihcation through the Majorana Fermions 
have been observed na. 

The purpose of this paper is to present a geometrical formulation based on the momentum 
representation of the coordinates in periodic solids. This method has been introduced in 
the context of the Spin-Hall effect [12] . In order to study the physical properties of periodic 
solids we use the Brilouine Zone {B.Z.). Due to the spin degrees of freedom the wave 
function is replaced by a spinor. When we compare spinors at different points in the B.Z. 
we observe that the spinors rotate in the B.Z.. For this reason we need to introduce a 
method which compare the spinors at different points in the B.Z.. The parallel transport 
[T3] method is best suited for such problems. In momentum space the coordinate x is 
given by the momentum derivative idk. Since the spinors depend on the momentum the 
derivative will be covariant and will introduce a spin connection. These comparison of 
spinors at different points in the B.Z. is done with the help of the spin connections (gauge 
helds) [l2]. For different cases we have either a multivalued spin connections, obstructions or 
degeneracies PI- The discrete symmetries: time reversal symmetry, parity inversion, mirror 
symmetry or charge conjugation acts as constraints which £x the topological properties. Due 
to the symmetries the eigenvectors satisfy certain constraint equations. The solutions of the 
constraint symmetry gives rise to specihc gauge symmetry for the spin connections in the 
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momentum space. This gauge symmetry is used to compute the electromagnetic response, 
determine the topological invariants . The electromagnetic (magnetoelectric) response for 
the T.I. is characterized by the second Chern number. For the Weyl Semimetal one of 
the symmetries, time reversal or inversion are broken and as result one has an even number 
of Weyl nodes connected by Fermi arcs. The Berry curvature for the Weyl nodes are 
at opposite momentum and gives rise to monopoles and antimonopoles. 

The interplay of Topological Insulators (T.J.), Weyl metals and Superconductivity 
gives rise to Majorana Fermions 

The plan of this review paper is as follows: In section II we introduce the method of 
covariant coordinates in the momentum space. We show that this approach gives rise to the 
concepts of spin connections and curvature. In section III we use the method to study the 
surface properties of a TJ. We show that the spin connections for = — 1 gives rise to 
anti-localization. Section IV is devoted to the Heisenberg equation of motion given in terms 
of the covariant coordinates . Section V is devoted to the computation of the Topological 
invariants for the T.I.. In section VI we compute the Topological invariants for Topological 
Crystals. Section VII is devoted to Weyl Semimetals. In section VIII we investigate a T.I. 
for a curved surface and show that the eigenvalues are effected by the curvature. Section 
IX is devoted to topological invariant for Superconductors. In section X we consider the 
p — wave superconductors which emerge on the TJ surface . Section XI is devoted to the 
Andreev crossed reflection and computation of the differential conductivity. In section XII 
we show that the Majorana fermions can be detected with a Squid in a persistent current 
experiment. In section XIII we present our conclusions. 

II -The method of parallel transport in momentnm space, spin connection and 
curvature -essential tools for Topological Insulators 

Topological Insulators and Superconductors are gaped fermionic systems which exhibit 
topological protected boundary for an arbitrary deformation, as long the discrete symmetry 
such as time reversal, particle hole and chiral symmetry are respected. . Due to the 
symmetries the Hamiltonian in the momentum space h{k) is invariant at the time reversal 
invariant point k = T (points in the Brillouin zone) which obey Th(F)T“^ = h(—F) (T is the 
time reversal operator which obeys = —1) or charge conjugation Ch{k)C~^ = —h*{—k) 


3 



(C*^ = —1). The product of the charge conjugation (particle-hole) with the time reversal 
allows to define the unitary chiral symmetry which holds in the entire B.Z.. As a result 
we have 10 symmetry classes. For the case that the inversion symmetry P and the time 
reversal T symmetry hold, it has been proposed [HI [58] that computing the eigenvalue for 
the ground state parity state, = SgnM(r) (M(r) is the gap at the 

momentum F) determine the topological invariant index Z 2 . This index is given by the 
product / = n.SgnM(F,) of the parity eigenvalues at the time reversal invariant point 
Fj in the B.Z. The index I determines the coefficient of the quantized electromagnetic 
response [ 231121 I- 

The spinors in the presence of the spin-orbit interaction vary from point to point. We 
will work with the gauged transformed Hamiltonian h{k) defined as = 

h{k,idk + A{k)). A{k) is the spin connection which emerges from the covariant derivatives. 

For translational invariant systems the Bloch spinor are given by \k Zi fjaik)) = \Ua{k)) 
(when orbitals are included \k Z (Ti^k) Z Ts{k)) = \Ua{k)) ). As a result we have the repre¬ 
sentation, 

\m} = ( 1 ) 

a 

When the coordinate acts on the spinor we obtain: 


x^Wk)) = ta„ji,(k)) = ( 2 ) 

a 

a is the index of the band and includes also the spin degree of freedom. Due to the spin 
connections the derivative is replaced by the covariant derivative: 


- tAy{k)c^^\u^(k)) 


A,{k) = A“J{k)\U„0)){U„(k)\, -iAl-^(k) = (U„(k)\d^\Uf(k)) = 

Xa = Xa + Aaik) 


dxU:^,{x){-tX^)Up4x) 


( 3 ) 


Xa is the coordinate in the momentum space, Xa = idko. which obeys [a;“, fc;,] = [X“, kh] = i5a,b- 
The covariant coordinate is defined by Va = Xa = Xa + Aaik). 

The covariant derivative introduces the concept of parallel transport, which shows 
how spinors are translated from point to point. The parallel transport condition is 
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{^l]{k)\Vai^ik)) = 0 which gives rise to the equations 


\^{k,P))=P 


.-*/-oo dk'-Aa{k') 


\Pik)), == P 


dk'-Aa{k') 


mm (4) 


Where P defines an ordered path. From the spin connection operator Aa{k) |T3] we construct 
the curvature operator- F{k). This is done with the help of the exterior derivative d. Acting 
with the operator d on the eigenvector \Ua{k)), id\Ua{k)) we obtain the representation of 
the spin connections Aa{k) = 'Yl,otidAa'^{k)\Ua{k)){Uii{k)\. The curvature operator Fapk) 
is obtained from the spin connections A^'^{k). 


F,.,Ak) = [w(i),xt(fc)] = 


a,l3 


A/A) = aAfHk) - d,AZ-Hk) + i[daAi-\k),d,x'/(k)] 


( 8 ) 


The curvature allows us to introduce the first and second chern character |13]. 

chi = ^Tr(F,,b(h)), ch2 = ^e“’'’‘=’''Tr(F,,fe(h)Fe,,(h)) (6) 

In Appendix A we will show how the curvature is derived from the eigenfunctions. 


Ill- Observing the topology using external sources or disorder potential- 
transport of Topological Insulators 


The Hamiltonian can be express in terms of the eigenvalues. To obtain information about 
the topology, we have to transport the spinor around the B.Z.. Alternatively we can couple 
the scalar potential ao(x) to the electronic density p{x) and the vector potential a(x) to the 
current density J{x). 

For spin half Hamiltonian with two orbitals \k®a =t, i = 1, 2) for the T.J. Hamiltonian 
Hq we find: 

H = Hq + J d'^x[p{x)ao{x) + J{x) ■ a(T)] = Hq + (7) 

The four component spinors for the Hamiltonian Hq\ 

p{x)= f dvmik), p{k) = j^mmusik)+Bi{-k)Vs{-k)] (s) 

s=l,2 
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Us{k) is the four component spinor for the particles and Vs{k) is the four component spinors 
for the antiparticles, s = 1, 2 represents the spin helicity. 


= j d^qY,[cl{k)V{q){Us{k)\UAk + q))CAk + ^ 

s,s' 

[ d'^kCl{^6s,snda + Ai^'^'\^CAk) 


s,s' 


t6s,s'da + Ai^'^'\k)=X^, 


X\X^ 


= FaAk) 


(9) 


On the surface of a T.I. (Topological Insulators) the T.R.S. (time reversal invariance) gives 
for the spin-orbit scattering, the integrated Fermi Surface curvature is tt which gives rise 
to anti — localization. This can be seen from the properties of the spin connections. The 
spin connections change sign for the transformation k —)■ —k. This change of signs occurs 
also for the Cooperon channel! As a result the conductivity increases. The two dimensional 
TI Hamiltonian with the scattering potential V{x) is given by: 


h{k, x) = -a 2 ki + aik 2 + V(x) 


( 10 ) 


For a hnite chemical potential fa > 0 we consider only the conduction band. We hnd in 
terms of the eigenvectors UAk) the representation: 

^Ak) = C{k)UAk) 

U„=Ak) = x(-fc) = xA) + tt 

( 11 ) 


We hnd for the conductions electron the scattering potential gives rise to the Hamiltonian: 


j^ext. _ 


(i^a;T''’(a;)H(T)T(T) = E d^k / dA\v{^C\k){UAk)\U„ik + ^)C{k + ^ 


/ d‘^qV{q)){-iqA / d^kC\k) ida +Aa(k) C{k), daAa(k) = da-x(k) 


( 12 ) 


Eq.(12) determines the scattering matrix S. We hnd that the scattering matrix S obeys: 


S{k -k) = A^TS{k -k) = -S{-k k) 


(13) 
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IV-Equation of motion in the B.Z. for non commuting coordinates -an exact 
formulation 


In this section we will use the Heisenberg equation of motion using the representation 
introduced in the previous chapter: Xa = idk^ , [x“, fcb] = [X“, k}] = Xa = Xa + Aaik). 
In addition, for a non uniform system we will have a real space curvature similar to the 
magnetic held. 

The methodology of a ’’curved” space induced by the spin orbit coupling in the B.Z. was 
described in terms of the spin connections = {Ua{k)\ida\Up{k)) and the curvature 

The Hamiltonian in the eigenvalue representation is given by: h{k) = '^xEx{k)\U\{k) >< 
U\{k)\ the presence of a scalar potential V(r) is replaced by V(i?) where R is the covariant 
coordinate. We hnd for the Heisenberg equation : 


h{k,R) = h{k)+y{R) 

IF - 


( 14 ) 


When the real space is curved, due to dislocations or magnetic helds the momentum operator 
is replaced by a covariant momentum = k°‘ — ^Log[g] where g = x/Delg^ is the 
metric tensor na. For this situation we introduce the real space momentum curvature 
fla,b = [K°',K^] which together with coordinates curvature Fa^ = [Xa,Xb] modihes the 
equations of motion. 


giFo. _i ^ 1 

= -^dnX{R) + 

at n n 


( 15 ) 


V-Topological invariants from response theory 

We will demonstrate how the topological invariant appears. In this chapter we will appre¬ 
ciate the power of the covariant coordinates and their commutations. We observe that the 
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topological invariants are given by the commutators of the covariant coordinates. 

We will consider a typical T.I. Hamiltonian for the materials Bi 2 Se^, Bi 2 Te^, Sb 2 Te^ . 
We introduce the tensor product |a) = \a =t, 4-) ® |t = 1,2) (cr stands for the spin and 
r stands for the orbitals). A four band model is obtained |62] which can be written in the 
chiral form |37] : h{k) = Ie(k) + 71^2 + 72^1 + 737^3 + loM{k). 

The first term affects only the eigenvalues and not the eigenvectors. The parameter rj 
satishes p << 1. The 7 matrices are given as a tensor product ■■li = cr* O (-l)r 2 , i = 1,2,3 
, 7o = / 0 Ti ,r 5 = J X r 3 . The mass (gap) M{k) obeys M{—k) = M{k) and has points in 
the Brillouin Zone where it vanishes. (On a lattice with the lattice constant a we dehne the 
Cartesian component of the momentum ki = TPe Hamiltonian h{k) is diagonalized 

using the four eigenvectors \Us^\k)), s = 1 , — 1 are the spin helicity operator and e = +, — 
represents the particles-antiparticles energies, E{k) = e{k) ± \Jk'^ + M‘^(k) with the mass 
M{k) which vanishes at k* and M(r) 7 ^ O.The Hamiltonian in the eigenvalue representation 
is given by. 


h{k) = 5 ^ [B(«)|f/i'-+i(b>(7"*>(^i + E(k)\ui'‘=*Hk)){u‘r*\k)\ 

s=l,-l ^ 


(16) 


The Green’s function operator in the a basis |a) = \a =t, i) 0 |t = 1 ,2) is given by: 
G{U, k) = <^(^1 fc)a,a'|«)(«'|- 

In the eigenvector basis the Green’s function takes the form: 

G(.. k) = y ^ (,,, 

oj — E(^k') id C(2 T £'(— k') T ia 


The transformation from the |ct:) basis to the eigenvector basis \ Ug{k)) replaces the coordinate 
fi = idl with the covariant coordinate Ri [12]. In the second quantized form the spinor 
operator T(r) is given by: 


T(f) = 


d'^k 

(27r)« 


Ak-r-y 


'^(k),'i!{k)= y \c,{k)uf+\k)+bi(-k)up''{-k)\,'i!{k) = ¥{kho 


s=l,-l 


(18) 


(It is important to stress that this representation is valid for momentum \k\ < \k*\, for the 
region \k\ > \k*\ we need to choose a different representation.)The coupling of the T.I. to 
the electromagnetic field d{r,t) and ao{r,t) is given by the action S'®’”*: 












3 3 


u=0 11=0 

^{k,u) = '^'<{k,u)'yo]Rf,=o = id^ 




(27r)'^ J 27r 


(19) 


We compute the partition function [37j for four space dimensions. We find that the 
effective action for the electromagnetic fields 0 is given by r[ao(t, f), d{t, f)]: 

2 =/nnn dCs{k,uj)dCl{k,u)dbs{-k,u)dbl{-k,u)R^^"+^^^^ = e^n-odRAtxl] (20) 

S k iO 

Using the totally antisymmetric tensor ecn,a 2 , 03 ,« 4 ,a 5 four space dimensions we find the 
electromagnetic response which is given in terms of the electric E and the magnetic field B. 
The coefficient C 2 is quantized, for values oi 6 = 7i,0 when the time reversal symmetry is 
preserved |H]. 

The electromagnetic response in 4 + 1 space time dimensions is 0II31I231E11I151I79]: 


C 2 = const. k)RoG-\k)G{uj, k)RiG-\uj, k)G{u, k)RkG-\uj, k)G{u, k)RiG-\u, 

( 21 ) 

In the presence of electron-electron interactions the Green’s function can have zero’s 
[541 [55], for such a case the system is not topological! If the renormalized Green’s func¬ 
tion has no zero’s the topology is preserved. The renormalized Green’s function 
Gr{u:, k) = E. 


tion Z{j,Ui,M) = 


^ (uj+ie—Es{k)) 


1 - 


\Us{k)){Us{k) \ is given in terms of the wave function renormaliza- 
- -1 

UJ = K 


du 


When the wave function renormalization Z is finite at a; = 0 we take the limit w —)■ 0 and 
obtain the second Ghern character [55], Z cancels: 


C 2 = constR^^^’‘^^Tr[G{k)R^G-\k)G{k)RjG-\k)G{k)RkG-\k)G{k)RiG-\k)] ( 22 ) 


The trace operator acts only on the occupied bands. The Green’s function in the eigenvector 
basis representation makes use of the covariant matrix coordinates i?*. The Ghern character 
is given as a matrix product. The commutator [Ri,Rj] = Fij{k) gives the curvature Fij{k) 
in terms of the spin connections \k). The second Ghern number is given by, 

G 2 = ^ / d^ke^'^''^'^Tr[R,{k)Fk,iik)] (23) 
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The integral in Eq.(23) is either zero (exact form) or non-zero (non exact form). Therefore, 

Tr[F,,,{k)Fk,i{k)] = Tr[F^] = d[K,] (24) 

where = Tr[AdA + \A^] is the Chern-Simons three form [13] . can be found with 
the help of the gauge symmetry imposed by = —1 [13]. The second Chern number in four 
dimensional space is given by: C 2 = 3 ^ which has a Z 2 winding 

number. Tr[F‘^] is dosed but not exact [T3|. This means that in some restricted regions of 
the Brillouin zone the integral J^ 4 Tr[Tr[F^] is given by a Chern-Simons contour integral 

m- 

Tr[F2] = = d{Tr[AdA + \a^]) (25) 

o 

This result does not hold in the entire B.Z.. We can identify two regions which are related 
by a transition function (a gauge transformation), a transformation matrix between states 
for the region k which belong to half of the positive sphere (5*^)+ and k which belong to 
second half of the negative sphere (5^)- . 

Next we compute the effect of the transformation fc —)■ — A: on the spin connections. Due to 
the time reversal symmetry we hnd the matrix which transforms between the two regions of 
k . The transformation is given by the PfafRan matrix B dehned in terms of the Kramers 
pair. We use the eigenvectors to compute the matrix Bs^s'{k) (T = —ia 2 K, K is the complex 
conjugation operator), 

B.Ai) = {Ui-H-k)\T\Ul,7\k)) ( 26 ) 

The relation between the spin connections A{k) and A{—k) is given by: 

A{-k) = B{k)A*{k)B\k) + iB{k)dk^B\k) (27) 

As a result the curvatures transform like : 

F„(-fc) = -B{-k)F:/k)B{-k) (28) 

Applying Stokes theorem for dlK^] in the two regions one obtains a boundary integral over 
the difference of the Chern-Simons terms (dehned for each region). The difference between 
the two Chern-Simons terms can be understood as a polarization P{q) difference between 
two regions. (The boundary term is a surface perpendicular to the fourth direction q. ) 

^2 = i = [ {[Tr[A+dA+ + ] - Tr[A.dA. + IaI]) 

JbZ JdS* 
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= ^ [ c^'Ki,fcTr[(5(fc)a,St(fc))(S(fc)a,i?t(fc))(5(fc)afci?t(fc))] = J2Nk*= 2P[q] 

^ k* 

(29) 

We have 2P[q] = N^* [k* are points in the Brillouin Zone where the Pfaffian matrix 

vanishes) . Due to the lattice periodicity, Bloch theorem allows to define the polarization 
F(^ as modulo an integer. We recover the result P = 0 for Ylk* = even and P = \ 
for Ylk* = odd [TT]. To conclude P represents the topological invariant for d = 3 space 
obtained from response theory. 

VI-Topological invariants for Topological Crystals 


The method used in the previous section is applicable to Topological Crystals where a 
mirror reflection invariant rj = + o' 2 )K) with the property 7f = —1 replaces the 

= — 1 invariant. 

Inspired by [58] the authors in ref. [59] proposed that SnTe has a mirror plane perpendicular 
to the [110] direction. A band inversion at the four L points in the Brillouin zone between 
SnTe and PbTe can be achieved for the mixed crystal Phi^^Sn^Te. The k ■ p model near 
an L point [59] gives the Hamiltonian: 

= (cr2 ® Ti)ki - {ai (g) Ti)k2 + (/ <g) T2) + M{k){I ® T3). (30) 


where M{k) is the inverted mass (has zeros in the Brillouin zone), a = ±1 corresponds to 
the states with total angular momentum J = 1 , |j = 1 , = ±1 > and r = 1, 2 corresponds 

to the p orbitals of the cation (Sn or Pb) and anion Te. 

The mirror invariant p with the property 77 ^ = —1 can be found such that the condition for 
the polarization is different from the one given by the time reversal invariant points. 

The reflection symmetry from the plane perpendicular to [1,1,0] is given by the transfor¬ 
mation [fci,/c 2 ,^ 3 ] [—^ 2 , —^ 1 , ^ 3 ]- The operator of reflection which acts on the states is 

given by a rotation of an angle tt around the axes [1,1,0] and is accompanied by an in¬ 
version trough the origin. A simple calculation shows that the mirror operator is given by 
M = -|- < 72 ) 0 /]. As a result the state |0(fc) > is transformed to \((>M{k') > and the 

Hamiltonian h^{k) obeys the symmetry: 


|0M(-fc2,-^ 1 ,^ 3 ) >= M|0(A)1, ^ 2 ,^ 3 ) >= + (^ 2 ) ® I]\(j)iki, k2, ks) > 
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M ^h^{ki,k2,k3)M = h^{-k2,-ki,k3); 


M = 



[(o-i + (J2) ® /]. 


(31) 

Next we include the anti unitary conjugation operator K and define the operator rj = MK = 
(;^)[((Ti + (T 2 ) ® I]K which obeys rf = —1 (this is similar to the time reversal operator 
T = —ia 2 K). We obtain the invariance transformation: 


?7 ^h^{ki,k2,ks)r] = {h^{-k2,-ki,k3)y. (32) 

At the invariance points k^ one obtains the conditions: 

uL!_UM _uM _uM\ _ uLf_uM _UM uM\ 

lb [ ) — tb y bv 2 1 1^'i )■ 

For two dimensions we have only two invariant points [0,0] and [tTjTt] and for three 
dimensions we have four invariant points [0,0,0], [0,0, tt], [tt, tt, 0], [7r,7r, tt]. 

At this stage we will follow the strategy proposed in section A. We identify the pairs of 
degenerate eigenvalues. For this case we expect to find a Kramers pair for the eigenfunctions 
|W=i(A;; M(/c)) > and | W=_i(/c; M(/c))). When the mass parameter M(k) has zeros we 
observe that the matrix fW,s'(fc) defined by Ws^s'(k) = (Vs(—k;M(—k))j7]jVs(k;M(k))} is a 
PfafRan which vanishes at some k , M{k) = 0. Due to the zero’s of the Pfaffian it is not 
possible to construct a single eigenfunction for the entire Brillouin zone. We observe that 
the eigenfunction at —k is related to the eigenfunction with k in the following way: 

|K=i(-fc;M(-fc) >= WU^^,^_yk)\Vs>=S-M{k) > 

\Vs=.i{-k-M{k) >= (fc)IK'=i(fc;M(fc) > 

(33) 

where M{k) = M{—k). For each pair of mirror invariant states at the points [0,0] and [tt, tt] 
(for d=2) and [0, 0, 0], [0, 0, tt], [tt, tt, 0], [tt, tt, tt] (for d=3). The PfafRan matrix W induces a 
transformation for the spin connections, resulting in the condition for the Chern-Simons field 
polarization 2P[q] = YlkM (with the sum restricted to the 77 invariant points). It results 
in a different condition for the polarization, since is different from the condition 

Ylk* for the time reversal case T = —ia 2 K. Therefore we can have a situation where 
the polarization is zero according to the time reversal symmetry and non-zero according to 
the mirror symmetry. 


VII- The Weyl Semimetal 
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In this section we will consider Hamiltonians of the form H = vpa ■ b{k) where b{k) is 
a vector-valued function b{k) has singular points kr where b{k) vanishes. Acording to the 
Nielson Ninomiya theorem the singular points kr are even and satishes ±kr- We dehne the 
Hamiltonian around these points: „ H{k = kr + k,p = —kr + p) = X]r=i n ' P- 

Around each singularity we have a monopole or an antimonopole. The sum of the monopoles 
and atimonopoles must be zero in a crystal. 

Next we will consider explicit examples: 

The Dirac metal {DM) has two Weyl nodes at a single singularity: H = -^vpa ■ k . When 
we dope the DM becomes Dirac semimetal {DSM). 

A Weyl Semimetal is characterized by an even number of nodes with an even of point 
singularity. The Hamiltonian has two opposite chiralities : H{±.kr + 'p) = ^vpa-p — p. where 
fi is the chemical potential. For each chirality we can compute the eigenspinors \Ua{k)) with 
the ground state energy a\k\. For each chirality a = ± we compute the spin connection, 
curvature and the divergent of the curvature : 

-tA:{k) ^ (uMdalUM), x: = Xa + A-{k) 

F:^,{k) = daAUk) - d,A:{k) 

Q;27r(5(/c CX ib 

( 34 ) 

For the monopole and antimonopole at ±kr- 

The time reversal symmetry and the inversion symmetry gives a restriction for these cur¬ 
vature and and the monopoles density. Either inversion symmetry or time reversal can be 
broken but not both: For broken inversion we have for the curvature, Fafi{k) = Fa^b{—k). 
For broken time reversal the curvature obeys, Fa^b{k) = —Fa^b{—k). 

On the surface, of a crystal the two Weyl nodes are connected by Fermi arcs, for this we 
have to demand that in the normal direction of the surface the wave function is zero. 

a) Realization of the Weyl Hamiltonian for broken inversion symmetry 

We introduce the Pauli matrix r for the two orbitals and a for the two spins: The model 
which breaks inversion is given by: 
h{k) = vpT^ ® a -k + Ti® IM{k) -f rs 0 Ig{k) 
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M{k) = M{—k) and g{k) = g{—k) are even functions. The spinor are given by: T(fc) = 

r 1 ^ 

We compute the eigenvalues and hnd the points in the B.Z. where E = 0. 

According to general theory the low energy Hamiltonian has at least two Weyl nodes H = 
±Vp^a-p. with the connections Aa{k-, kr), Aa{k-, —kr)) and monopole charge : ea,b,cdaFb,c{k) = 
±27r5("fc — (±fcr) 


b) Realization of the Weyl Hamiltonian for broken time reversal symmetry 


Superconductivity with broken time reversal symmetry has been considered in the literature 
by 0E]. Here we propose a continuum two band model with broken time reversal symmetry 
(T.R.S.) given by: 

h^^\k) = vf{ti (g) 0 - 3 )(/c^ - M) + vf{t 2 ® o- 3 )k 2 + vf{t^ ® o- 3)^3 + {I ® 0 - 3 ) 8 , M > 0 (35) 

The model has monopoles since the inversion symmetry is not broken. This model assumes 
that in a magnetic held B the two bands are polarized and the coupling is linear in the 
momentum space (spin orbit interactions and uniaxial stress can realize such a model ). 
The Pauli matrix r acts on the orbital space and the Pauli matrix a acts on the spin . We 
will project the Hamiltonian on the polarized state \U^{k)) and hnd: 

j^{Project)0^ = VpTiikl - M) + + VpTsks - I\B\ (36) 


In order to study superconductivity we consider an attractive short range pairing interactions 

^paring. 


=-H(|d|) / 


'iW(x+ ^)>i>(£+ i) + <!-*(£+ 


where \k(a;) = pi{x),ijj 2 {x) are the spinor buid from the orbitals 1 and 2. The paring order 
parameter must be odd to ensure the inversion symmetry. We use the Hubbard Stratonovici 
transformation and rewrite Eq.32) in momentum space for a hxed pairing distance d. 
■|A(d)p 


H, = 


d^k 


+ ['ilj\{k)p\[—k) + 'ilj\{k)'ilj! 2 {—k) jiA(d) sin(fc ■ d) + h.c. 


.2V(d)| 

We linearize the spinor around the monopole and anti-monopoles, g = 


(38) 


± \/M, k2 = 0, ^3 = 


0 


As a result we hnd the chiral Hamiltonian: 

hiik = [\/M -f Pi,P 2 ,P 3 ) = 2^/MtiPi - T2P2 - T3P3 
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hnUi ^ [-VM + Pi,P2,P3^ ^ -2'/Mti'Pi - T 2 P 2 - 

( 39 ) 

Next we expand the spinor in terms of the chiral fermions Cl{p) and Culp) determined by 
hiik) and hR{k). As a result the two component spinor is represented in terms of the 
chiral fermions . "^{p) = Cl{p)\V{p) + Ci?(p)|f^(p) • |ld(p) is a two component eigenvector 
which represents the two orbitals with left chirality, similar \U{p) has two components for 
the two orbitals with the right chirality. Using the eigenvectors |U(p) ,\U{p) we project the 
pairing held in terms of the two chiralities : 

i’\(M(-P) + 4(p)4(-p) ~ Cl(p,C'„(-^([l\V(pt )'+ (2\V(P)'([2\U(-P)') 

+U(p)cl(-p)({iV(p))*{i|r(-pT + {2|(7(p)-{2V(-p)-) 

= Ci(p)C'^(-p)F(p^ - Ci(p)Ci(-p)F(-p)-,F(p) = {lV(p-))-{l|(7(-pT + {2\V(p)-(2\U(-p)- 

( 40 ) 

Where (1| and (2| are the two orbitals. The Hamiltonian in the chiral form with the positive 
chemical potential p > 0 is given by: 

H = I d^p[cimf^-k^)CL{p)+cumm-f^)CR^^ 

+ [cl{p)Cl{-p)F{:^ - C'i(p)C'J:(-:p)F(-:^)iA(d) Sin(p- d) + h.c. + ^7^ 

(41) 


We construct the BDG Hamiltonian in terms of the spinors 

^. The pairing order parameter is represented by 

A(p) = iA((i)sin(p ■ d)F{p) which obeys A(p) = —A(—p) , A(p() = |A(p)|e*'^*^^ with 


ip{-^ = (p(p) + TT. 

From Eq.(27) in the previous section we have: C\{p){^i — 
ia^{p)^CR{p). 

We perform a gauge transformation where P represents the path order operator Cl{p) = 


P Cl{p) , Cr{p) = P Cr{p) The phase of the order parameter 

transform as ip{p) —)■ (p{p)—A{p). As a result the gauge held transforms as {d^{p)—d^{—p) —)■ 


d^{p) — d^{—p) — dpA{p) . Therefore the pairing phase obeys the equation : 

dp(p{p) - (d^{p) -d^{-p)^ where (d^{p) - d^{-p)^ = -p). (p; 

obeys: f —p) ■ dS = 47rg with the monopole charge q = 1- 
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The order parameter is given by: A(p) = |A(p)|e*‘^® , with the phase (p{p) to be determined 
from the solntion dip{p) = —p). Due to the multivalued of the pairing potential 

the phase can not be determined! 

The monopole charge allows for an expansion in terms of a set of complete eigenfunc¬ 
tion given by 0(p)) [ZS|- For a monopole it is not possible to choose a sin¬ 

gle vector potential in the entire space, one choses a vector potential in one sector. 
The complementary section is obtained by the gauge transformation As a result 

also the order parameter must be defined as a section. The Monopole Harmon¬ 
ics Theorem allows to expand any function which obeys the section condition HB]: 
|A(p)|e*‘^® = m where the complete set of eigenfunctions 

yq,i,m{0{p),(l)ip)) are given by, 


Yg^l^rn{ 0 {p),((){p)) = “ COs{9) 


(-q-ra) 


1 -f cos{9) 


(q-m) 


P, 


(—g—m),(g+m) 


l-\-m 


cos(9))e‘‘’"+''>* 

( 42 ) 


From the formula for Yg^i^rn{9{p), (j^ip)) we find: 

Yg=i^i=o^m=o[9{p},(j){p)] « sin(6')e*'^,Fq=i,/=o,m=i[6'(p),0(^] oc cos^(6')e*'^, 
Yg=i^i=o^rn=-i[9{p),(j){p)] oc sin^(6')e*^'^ 


(43) 


This shows that the order parameter A(p) = |A(p()|e*‘^® = 

m (^q=i,i,m{\i^)Yq=i^i^m[9{p), 0(p)] is different from the p — wave order parameter 

c) The effect of a scattering potential in a Wey material with broken time 
reversal symmetry 


Here we observe the advantage of our procedure given in Eqs.(9,12) for a scattering 
potential V{x). 

f dW ^))(-0 [ d^P^^^ip)[lda + A^ip;K)\^R{p) 

j ida + Aa{p]-kr) 

CL 

( 44 ) 


Where A(p; k^) is the monopole vector potential and A(p; —kr) is the antimonopole vector 
potential. 
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For a monopole it is not possible to choose a single vector potential in the entire space, one 
choses a vector potential in one sector. The complementary section is obtained by the gauge 
transformation [TS] where the charge of the monopole antimonopole is g = ±|. As a 
result also the order parameter must be dehned as a section. A theorem named Monopole 
Harmonics Theorem says that any function which obeys the section condition can be 
expanded in the complete set of Monopole harmonics. The complete set is given by : 


^=±TZ,m 


{0{k), (j){k)) = Mq^i^rn ( 1 “ COs{9) 


{-q-m) 


1 + cos{9) 


{q-m) 


P, 




l+m 


cos(6*))e 


i(m+q)4i 


(45) 


Pi+m '^^’^‘^~^'^\cos{9)) are the Jacoby polynomials [78] 


As a result we have the expansion for the spinors and 






L,m 


l,m 


(46) 


This representation can be used for computing the scattering effect in the presence of dis¬ 
order. The advantage being that using this expansion we do not need to work with the 
singular vector potentials e A(p; ±kr) 


Vlll-Surface states for an arbitrary crystal-face boundary 


Photoemmision, photoconductivity, optical conductivity and scanning tunneling mi¬ 
croscopy are sensitive to the nature of the surface states. The topology of the surface 
states is affected by the physical boundary. For an arbitrary surface we will use a curved 
coordinate basis. When the coordinate system rotates from point to point we can use a 
non-coordinate basis, da, a = 1, 2, 3 introduced by Cartan nsum which is related to the 
Cartesian system by the vector , p = x,y, z. For example a point on a surface is given by 
the cartesian coordinates r = [x{ui,U 2 ),y{ui,U 2 ), z{u\,U 2 )] where {ua=i,Uh= 2 ) are the coor¬ 
dinate on the surface dehned by the normal Ej^ to the surface dehned by {ui,U 2 ). The two 
sets of coordinates are related by the matrix = |^, y = x,y, z and a = 1, 2. The normal 

dr ^ dr 

to the surface is given by E^ = This set of transformations allows to replace 

11 ^ 11 


17 



a^'da by the covariant derivative Va = 5a + dc] (which depends on the connection 

dehned below) ; cr“Va = 5^ + |r“’^e))[(Ja, o';,] , and e); are the transformation and 

the inverse transformation matrix and is the connection one form matrix (see 
page 285). 

The Weyl Hamiltonian in cylindrical coordinates has been consideredin the literature 
[701 [731 [H]. Here we will modify the method in order to deal with the arbitrary crystal face 
boundary. 

a) The Weyl Hamiltonian in Cartesian coordinates is = {—i)[axdy — ayd^]- 
We put the Hamiltonian on a cylinder and take the axes x = y = x‘^ = rsin[0] and 
z = x^ = rcos[0]; r = x,rsin[0],rcos[(/)] . We will study this problem using the non¬ 
coordinate basis [13]: Ua=i = x ,Ua =2 = 0 and r is the coordinate in the normal direction 
Eisf] the derivatives are da=i = da =2 = = dr and the differentials one form 

are given by: 9^ = dx] 6^ = rd<p-, 9^ = dr. The coordinate basis is not hxed, therefore 
connections r“ will be generated. The connections r“ j, are obtained from the Cartan’s 
equations. We have for the Torsion T“ the equation : A9^ = T“, a = 1, 2, 3. The 

connection is expanded in terms of the differential T^^9‘^ one form using the matrix e(): 

= 5“ ;E^ = daV = daX, dar sin[0], daV cos[(/)] 


= 


1 , 0,0 


■,E^ = 


0 , r cos[(/)], —r sin[ 


■,E^ = 


0 , sin[(;/)], cos[ 0 ]. 


The transformation will determine the form of the Hamiltonian in a rotated basis : and 


b,a 


, .a,b _ pa,6„c . , ,a,b _ , , 

From the torsion condition Torsion = T“ = 0 we hnd: d9°‘ + A 6 *^ = T“ = 0, ( A is the 
wedge product [T3|) a = 1 , 2 , 3. We obtain the equation: ;r 2 3 = —r |2 = 

As a result the Weyl Hamiltonian in the cylindrical basis is given by: 


a=l,2,3 


a=l,2,3 


( 47 ) 

The eigenvalue equation H^y’"‘ip{x,(j)) = E'il>{x,(j)) has real solutions for boundary con¬ 
ditions "0(^)0 + 27r) = —-^(x,0). We hnd: E = '0=‘=(x,0) = 

l^^-i[±x{k.,i)-a{k.,i)] ^ ,n = 0 ,± 1 ,± 2 ,..; I = 0 ,± 1 ,± 2 ,.. 


^ikj;X+{l+k)(l, 


;tan[x{K,l)] = ; tan[a{kr,l)] = 


2ryfc2 + 


i(i+l) ■ 


b) The Hamiltonian for the three dimensional T.J. with a cylindrical boundary 

We will consider a three dimensional T.J. with the mass dependent gap is in the z 
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direction. Such a T.I. has a surface boundary perpendicular to the axes z with crystal-face 
[x, y] plane localized ai z = L. To simplify the discussion, we consider a situation where 
the crystal-face is cylindrical with the cylindrical axis in the x direction and length > L. 
As a result, any point on the surface of the cylinder is given by the set of coordinates 
r = [x, Lsin((;/)), Lcos(0)]. 

The four-band model for the three dimensional T.J., Bi 2 Sez [m of size Lx X L X L is 
given by, H = + where r are the Pauli matrix for the orbitals and a describes 

the spin; 


H± = T 3 {-mo{z) - m 2 dl) -h 1 x 28 ^, H\\ = {-i)Ti{a 2 dy - aid^) - m\\T 3 {dl -F dj) (48) 


The mass is a function of z, niQ^z) = \mo\Q[—z + L] +{—M)Q[z — L], where Q[z] is the step 
function. For z > L, we consider M —)■ cx), therefore we obtain a zero mode on the surface. 
For a crystal-face which is cylindrical, any point on the cylinder makes an angle of 0 with 
the z axis. 

We consider a situation where the surface perpendicular to axis 2 ; (the direction of the mass 
gap is given by is is a surface of a cylinder. Any point on the surface can be viewed as 
rotated by an angle (j) (the new axis z' makes an angle (j) with the original axis ^ for which 
the mass gap has been introduced ) The axes z' becomes the radial direction on a cylinder. 
The transformed Hamiltonian iF|| is expressed in terms of the covariant derivative Va- H'_^ 
and Hu are given by : 


F/'_l = r3[-mo(-T^) - m2{ -t-t) cos^( 0 )d^] [1x2 cos( 0 ) - ixiaxsm{(j))dr 

cos(0) cos(0) 

-^11 = (-^)[ho-xCOs( 0)(-V2) + Ticry(Vi)] = ixi[ax 


cos(0) 




+ (Tydx + — J 


(49) 


The Hamiltonian has zero mode solutions H(r), = 0. The zero mode spinor 

H(f^ = r]r,ahl{f) is given by a product of the scalar function f2(r) and the spinor rjr,a ■ The 
scalar function is localized at r = L and vanishes for r —)■ cx). Using the eigenvalues of 
and Ti we hnd: 

hr=l,,T;cos(<«>0 = ® hr=-l,<T;co.(<A)<0 = ® ;7^[lw] ; ^ = ±1 

We dehne the vectors |1 >= ^[r = 1,(T = 1], | — 1 >= -^[x = l,a = —1] which allow to 
introduce the rotated Pauli matrices : 
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^3 = i[|l)(l|-|-l)(-l|], = |[|1)(-1| + |-1)(+1|] and 52 = |H|1)(-1H-^|-1)(+1|] 

We will use the eigenstates ?7r=i,cr;cos((?i)>o to compute the projections < Tia^ > |t=i < 'TicJy > 
|t-=i for the Hamiltonian ify = (—i)[ri(T 2 ; cos((;/))(—V 2 ) + i)] (see Eq.(48)) 

We hnd: 

< Tia^ > 1^=1 = 53 Cos( 0 ) ;< nay > |^=i = (- 52 cos( 0 ) + 5isin(0)) 

As a result the Hamiltonian hTy depends on the angle (p : 

jj^cos(^) ^ - *^^'^^^*^^ 53 ( 9 <jj,+cos( 0 )( 52 - 5 i sin( 0 ))( 9 a.-^(cos( 0 ) 5 i +52 sin( 0 ))] = 

(50) 

We observed that the rotation of the crystal-face for a cylinder is different from the surface 
state of a cylinder. In particular We observe a change of sign of the Hamiltonian for a large 
angles, 

To conclude, the eigenvectors determines the surface properties, such as spin texture, and 
surface currents . We believe that this results can be observed by photoemmision or hoto- 
conductivity. 

IX-Topological invariant for Superconductors 


For superconductors we can use the invariance in the entire B.Z. obtained by combining 
the time reversal invariance and the particle-hole symmetry. As a result one can find a 
unitary matrix T, T^ = 1 which anti commutes with the Superconductor Hamiltonian. One 
can show that the Hamiltonian can be brought to the form ( due to the Superconducting 
gap we can use a flat Hamiltonian): 


Q{k) 



( 51 ) 


From the relation c/i 2 = d[K^{A, F)] we identify the winding number z /3 = d[K^{A, F)]. 


z/3= / (F3[5l+,F+]-F3[5l_,F_]) 

J a s'^ 

(52) 

In one dimensions we have O = ^ /b ^ dkTr[q{k)dkq~^{k)]. 
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For two dimensions one identifies the Z 2 index with the Pfaffian matrix (the Pfaffian at the 
time reversal invariant point is equal to the matrix elements of the Hamiltonian components 

g(fc) [2SlEZllSn|): 


^=n 


pf 




(iet[g[f]] 


(53) 


See additional discussions in Appendix B. 


X-Superconductivity on a T.I. surface 


. We consider a situation where the chemical potential fi > 0 obeys /i > A (the pairing 
held). Following [51] we have for the the projected Hamiltonian, 


H = 


^(u|fc| — ij)c\k)c{k) + (Ac'*'(fc)e*^^^^c'*'(fc)c^(—fc) + h.c.] 


(54) 


where e*®*-*^^ = por the conduction electron the surface eigenvectors are : \u~^{k)) = 


V 2 


1, Aj gives rise in the momentum space to a vortex, A(A) given by F{-^) = 


' fco ' 


fco ' 




1 — e '“ 0 . The Hamiltonian in the long wave limit is obtained by replacing the kinetic energy 
hv{\k\ — kp) ~ ~ ^V\- The term ^^^^^-^(^) is replaced by which is a spatial 

derivatives in real space. The gauge transformation in the presence of a superconducting 
vortex is: A(r) —t ; A*(r) —t C{f,t) —)■ C^{f,t) —)■ 

In the polar representation the held C{r) = C{r^(f)) and C'l(F) transform in 
the presence of the vortex phase 6 {f), (^(r, 0 + 27r) = —C{r, 0) ,C'^(r, 0 + 27r) = —Cl(r, 0). 
Consequently for a 27r vortex the helds C{F) -iC^r) are double-valued. 

The chemical potential pe//. in the presence of a superconducting vortex at r —)■ 0 changes 
the sign. 


H = (Tr 


hv 


|A| 


CHr,(j))[—{-d‘^ - —dl)-fieff.{r)]C{r,(j)) - ^C(r,0)eA(d, + -d^)C'(r,0) 


2 ,kp 
i 


i 

—( 
r 


+ ^-^C\r,(t))e A(0^ - l0^)(0(r, 0) + - [ dt [ dV|A(f,t) 
2fco r ^ 9 J J 


(55) 


Folowing EH EE] we observe that this Hamiltonian has a zero mode Majorana solution 


H(r,0),H(r,0)J : [t/(r, 0), H(r, 0) 
the normalization condition, J < 00 . 


V* 


/—I 


iT 


The function obeys 

Wr ■» /'T 
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Next we consider few vortices, we limit ourself to iV = 1 (odd vortices) localized 
at the space time coordinates [Ri{t)], ....[R 2 nit)] . At these points the amplitude pair¬ 
ing held is |A(r, t)|. We consider a situation where the 2n vortices are far from each 
other. Each vortex corresponds to one Majorana Fermion [2Tj7j,i = l,..2n. A{f,t,) = 
|A(r, t, Ov{fR) = Ri{t), ....R 2 nit). The phase Oy{f,t) sat- 

ishes : Jo{f,t) = ^[8182 - d 2 di] 9 y{f,t) = Zi5[f- Ri{t)], 

J{f, t) = ^i[dtRi{t)]S[f— Ri{t)], Zi = ±1. In order to evaluate the effect of the Majorana 
vortices we introduce the Majorana density, 


n(f,t) = = ^ - R 2 a-i{t)) + 6{f - R 2 a{t)] 

a=l 


( 56 ) 


The Majorana current density is given by. 


n{f,t) = ^ ^dtR 2 a-i{t) 6 {f-R 2 a-i{t))+dtR 2 a{t) 6 {r-R 2 a{,t)) 

a=l 


( 57 ) 


The Berry phase involves the time derivative of the Majorana operators = 

^{R 2 a-i{t)R) and %a{t) = 7 {R 2 a{t)R)- We compute the time derivative 
and ^^{R 2 a{t),t) which are defined with respect the moving coordinates R 2 a-iit), R 2 aif) 
(the position of the vortices). We construct the dual action introduced in ref. [77j : 


S = I dt I d^rj{f,t) in{f,t)dt + in{f,t) ■ dr + ^e>^’''’^dubx{r,t)[-eAx{r,t) + ^al{r,t)] 




(58) 


In order to evaluate the effect of the motion of the Majorana vortices we replace the 2n 
Majorana fermions by n spinless fermions. The field acts as a constraint which 

imposes the relation = i^e^^u^xdvd^9y{fR), The vortices are 

described by |a)[(r, t). In a regular superconductor the vortices are screened and the ground 
state response of the external magnetic field gives rise to the Meisner effect. 

XI-Transport trough Majorana Fermions-Resonant Andreev crossed reflec¬ 
tion. 

a)-The Majorana modes 


22 








A direct evidence for a topological superconductor can be observed from transport. At¬ 
taching two metallic leads to a superconductor which support Majorana fermions one can 
observe the resonant Andreev reflection [22] and the Crossed Majorana fermion Andreev 


transmission . Superconductors have an Andreev reflection oc (£)^ which is smaller in com¬ 
parison to the resonant Majorana Andreev reflection which is of the order of 1. Most of the 
work done was for reflection ( an incoming electron is reflected as a hole and a charge of 2 e is 
propagating in the Superconductor) was investigated with the help of the S matrix [TTl l28l - 
122 ] . One of the difficulties for investigating the resonant Andreev crossed reflection (an 
incoming electron at lead 1 is reflected as an outgoing hole in lead 2 ) is due to the fact that 
the Majorana modes have a finite energy. The S matrix is obtained by imposing the conti¬ 
nuity equation an unitarity. The Majorana modes makes it difficult to compute the S matrix 


and to represent it as a Dyson series S = T 




where = j d'^rh{r,t) is the 


perturbing Hamiltonian. The knowledge of the S matrix given in terms of the Dyson equa¬ 
tion has the advantage of allowing the use of the ’’mashinery” of the Renormalization Group 
[271 ESI El]. The Majorana fermions have no coherent states,therefore the construction of the 
path integral looks impossible. This problem can be solved by halving the degrees of freedom 
and introducing spinless fermions which . For an even number of Majorana fermions 720-1 
and 72 a a = 1,2, 3,4... with the eigenfunctions F 2 a-i(F)) F 2 a{f) we have the representation 


for the zero mode fermions: Coir) = Co(r) = 


F2a{^ = 


1 pU2a 


—i 


hajr) 


^ and F 2 a_i(r) = 


72a-1^2a-l(h*) + l 2 aF 2 a{r) 
T 


, where 


1 p\4>2a-l 


1_g 2 *i^2a—1 

—i 


are the 
Vr 


two component spinors localized at the positions r = i? 2 a-i and r = i? 2 a- We introduce the 
spinless fermion operators Cl, Ca , 0 , = l,2,3...n and construct the Fermionic path inte¬ 
gral. The transformation between the two representation is given by: 720-1 — ^ Co + Ca 
’ F 2 a = ^ Co ~ Ca ; « = l,2,3...n Tile low energy Hamiltonian takes the form : 
g(Majorana) ^ J2a=i I Cl(*^t)Ca “ ^oClCa , ^o is the Overlapping energy for the Majo¬ 
rana Fermions ?eo72a-i72o- For an odd number of Majorana Fermions ( we will have for 
the 2n -|- 1 Majorana an unpaired Fermionic) ,we choose for 7^+1 = + (n+i 


or 


7n+l = 


iV 2 


Cn +1 ~ Cn +1 • The coupling of the Fermions to the two leads is given by 


Ft = tY^ f dy dl{x ==-^,y)+dl{x = = ^,y) + h.c. (59) 


a=t7 ' 


We consider two metallic leads at x = — F (first lead ) and x = ^ (second lead ). Keeping 
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only the Majorana Fermions we replace with the zero mode two component spinor 

(?o(r). 


b) The differential conductivity for the Andreev crossed reflection-A pair of 
two Majorana vortices 


Contrary to the single Majorana half vortex, the presence of two half vortices induces a 
current trough the superconductor. We will compute the Crossed Andreev reflection 
which is a process in which an incoming electron from lead 1, is turned into an outgoing hole 
in lead 2. In this case a single electron at each lead is tunneling into the superconductor to 
form a Cooper pair. For simplicity we consider two such vortices localized at the edges of 
the superconductor with the tunneling Hamiltonian and Majorana overlapping; 


Ht 




^(Majorana) ^ 


di{x 


L 

2 ’ 



7i + 




d 2 {x 


L 

2 ’ 



72 


(60) 


The overlapping of the two Majorana Fermions is given by ieo 7 i 72 which are localized at 
different position and are not-orthogonal. We replace the Majorana operators with their 
fermionic representation 7 i = C^ + C ^^*^72 = ^ C^~C • The Id leads are (ii(a; =—^) 

and d 2 {x = ^). We integrate the Majorana fermions and obtain the S matrix for scattering 
between the leads: 


S = 


dt 


dr 




J —oo J 0 ^ 

V^{t)V{t — T) = — di(t)e*^ +id2(t)e“*^ 

^ — d\{t — r)e“*^ + di{t — r)e*^ + idl{t — r)e“®^ — 


- id2(t)e"'^^ j ■ 
id2{t — r)e*^ j 


(61) 

We expand the fermion operators in the leads using the right and left movers. di{t){x = 
—1,0) = + d 2 {t){x = |,0) = R 2 {t)e^^^^ + We apply 

on the left lead a voltage H/2 and on the right lead a voltage —V/2. As a result we obtain 
for each leads two Green’s functions: 


Gj’^(E,a;) = 


e{E-f' 


{E-f)+t0 


+ 


u 


UJ 


d(-E + f) 


fO 
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Gl'\E,u) 

Gl'^{E,u) 

Gl'\E,u) 


9{-E+f) 
u + {E-f)-iO 
e{E+f) 
u-{E+f)+zO 
Q{-E-f) 
u + {E + f)-iO 


e{E-f) 
uj + {E-f)+iO 
Q{-E-f) 

^ u-{E+f)-iO 
GjE + f) 
uj + {E + f)+tO 


(62) 


Q{x) is the step function (zero for a; < 0 and one for a; > 0) The current in the leads is 
given by: J(a; = -L; h) = ev{N^E _ jyiT) _ _ -ev{N^’^ - In order 

to compute the current we will compute the Green’s functions as function of the tunneling 
parameter and the voltage applied on the two leads: G^’'^(i?, a;; ^), G^’^{E,uj; (left 
lead) and G‘^'^{E,u-, —^), G‘^'^{E,u] —^) (right lead). The self energies for each lead are 
a function of the frequency ;S^’^(a;), S^’^(a;) and S^’^(a;), S^’^(a;). 

We find, to order g‘^, the self energies: 


S^’^(a;) 

S^'^(cn) 

T{u,uo) 


/ 1 -U - _2_ 4 

-2T{ijj,uJo)—Ln[ -^ ) + f2T(a;,a;o)— sgn(a;) 

V Vi / V 


-2T{u, Uo)—Ln 

V 


■1 + 

2 

A 

T - 


A 

1 + 

-^ . 

A 

1 - 

W-^ , 


i2T{uj, cjo) —sgn(cj) 

V 


Tq + (cn + Wo)^ 


(63) 


Where A is the band with, Tq is a damping factor which is and Hujq = eo is the Majorana 
energy. The imaginary part of the self energy obeys and 

the real part of the self energy obeys 3f?S^’^(a;) = 3f?S^’^(a;) = S^(a;). The Green’s 
funtions are given in terms of the self energies: 

= (u, - (E - 

G^’^(E, w; —) - (w + (i?-^ 

(64) 


The real part of the self energy is used to compute the wave function renormalization function 

Z. 



a ;=0 


z~i 
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( 65 ) 


z-' = 


1 + ? 


A Vl - ^ 

2A 


r = 


4:T{u = 0,Uo)g^ 




t.(r§ + a.§) 


f is the microscopic width dehned at the short lenght scale I = log =0 given in Eq.(64). 
Due to the nonlinearity of the effective action Eq.(60) the width E scales with the voltage 


I = logi \ ) . Using scaling equation 


for the coupling constant g^ we find 


that the the width satishes the equation: 


^ = —const.T^ with I = log 


eV 


We obtain the differential conductance for the Crossed Andreev reflection 


dl(V) 

dV 


= ^ de ^ 


dI{V) 
dV 

^ ^ / de^ 


e 
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-AZ 


'-AZ 


— -:: - { If D (D -p — Z) -p fpD — —Z 

27Ti{Q-ey + {rzy\dnV^-^-^ 2 ^ 2 ' 
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f(U) 


+ 


f(U) 


27rLr^_e)2 + (r(U))2 {f + ey + {T{V)y 


^ WcTan[^^(l + ^)] + ArcTan[^^ 


h TT 

r(/(u)) ^ r(u) = 




r(/ = 0 ) 


(1 + const.r{l = 0)Log[^]) 


,r(/ = o)^f 


( 66 ) 


We And that for a pair of vortices the Andreev crossed reflection obeys 
(see fig.l) and for a single vortex —)■ 0. 


v^o 


h 


XII-Detection of the Majorana Fermions for two metallic rings pierced by a 
magnetic flux which by p — wave 


We consider a situation where two metallic rings are attached to the two ends of the 
p — wave (which has two zero modes at the of the wire, 71 and 72 ) (see |22l SI])- 
We consider the special case where L = Nl {L is the wire length and I length of each ring ) 
The flux in ring one is /i and in ring two is f 2 - Using periodic boundary conditions Ci{—1) = 
Ci(0) and C2{—1) = (^ 2 ( 0 ) we perform a gauge transformation Ci{x) = e''^^^Ci{x), C 2 {x) = 
e^^f^C 2 {x) which introduce twisted boundary conditions. We And the spectrum of the 
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FIG. 1: the differential conductivity for the Andreev crossed reflection 


uncoupled rings as a function of the momentum n = 0, ±1, ±2,.... 

Ho = '^ ^C!(n)(n + /d^Cdn) + |c'i(n)(n + 


(67) 


{{Maorana _ jg Hamiltouian of the p — wave wire restricted only to the Majorana 
modes. We express the Majorana fermions in terms of the fermion helds ( and 


7i = 


72 1 


7 + C 


, 72 = T 


72L 


7-C 


( 68 ) 


The coupling between the wire and the two rings takes the form : 


^ E [(7'(«)(c + c') + (c + ctCiW)] + (-o4 E [4(")(c - C) + (c - c*)<? 2 (n)) 

^ n ^ n 

(69) 


We integrate the rings degree of freedom and obtain an effective Majorana impurity Hamil¬ 
tonian. 


H,n = c'^('^) V(u;/i,/2) C(u);C’'^(u) = C’M.C. 


( 70 ) 


M{u) is a 4 X 4 matrix which depends on A and 6 (the function obtained by integrating the 
ring degrees of freedom) . A and S are given by: 


A = f (Ad) + Ad)), <5 = f(Ad) + Ad)), Ad) = 


uj—Ei {n)-\-ix 


Ad) = 


E 7 


u—E2{n)-\-ix 

{Lj-E2{n))^-\-x‘^ 


,x —)■ 0. 
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The matrix M{uj) is given by: 

Mi^i{u}) = cu — eo — A, = 0,ikfi^4(a;) = —5; M2^i{u:) = 0, ikf 2 , 2 (^) = 

— (a; + eo)+A, M 2 ^ 3 {uj) = 5, M 2 , 4 ( 0 ;) = 0; M 3 , 1 ( 0 ;) = 0 ,M 3 , 2 (ci;) = 5, = —(a;+eo) + A, 

= 0; M4,i(c(;) = 6,M4^2{^) = = 0, = (cj + eo) — A; Eiiji) = 

f + hf - Ai; E2{n) = f (n + f2f - 

We integrate the Majorana Fermion and obtain the exact partition function: 


Z(/i, h) = z{g = 0; /., h) ■ det\M(ur, /i, /j)] 


(71) 


Z{g = 0) is the partition function for two uncoupled rings. The effect of the coupling is 
controlled by the Majorana contribution det[M{u:] /i, / 2 )]- Therefore the current is given by 


/l(/l,/2) = 


dLog(Z(/ij2)) 

dfi 


; h{fl,f2) — 


dLog(Z(/i,/2)) 

dh 


Due to the multiplicative form of 


the partition function the current is a sum of two parts = 0) (i=l,2) and a second 

part is determined by the matrix M(a;;/i,/ 2 ) and is given by 5/i(/i,/2) , i = 1,2 . The 
current in each ring is given by, 


h{fi,f2) — A(/i; S' — 0) + Mi(/i, /2; M) 
/2(/i,/2) = I2{f2;g = 0) + M2(/i,/2;M) 


(72) 


We investigate the case of equal fluxes, /i = /2 = /• Due to the fact that L = Nl the 
hoping matrix elements are real. In particular the Majorana energy cq couple like a regular 
impurity to a set of states determined by the two rings . Effectively the integration of the 
electrons in the rings renormalizes the energy cq to eeffX^o, /) = ^0 + ^(/) where E(/) is the 
shift in energy caused by the energy in the ring E{n, /). 

This result show that when the Majorana energy goes to zero the only contribution to 
the persistent current comes from the perfect wire in hgure 2 .In the presence of a Majorana 
Fermion we have in addition to the persistent current from the perfect wire the contributions 
^e//.(eo; /) = ^o + ^(/) determined by the matrix Mioj) given in £g.3 for the Majorana energy 
Co = 0.01 and in £g.4 for the Majorana energy cq = 0.1. The persistent current is measured 
by scanning with the SQUID |12] which measures the change in the magnetization (The 
magnetization is proportional to the persistent current). One subtract from the value of the 
persistent current, the current from a single ring , and hnd the dependence of the persistent 
current (the correlation part given by 5/j(/i,/2) , A = 1,2) on the Majorana energy. From 








I[fg=0| 



FIG. 2: The current for a single ring, the coupling constant to the wire is zero. The current was 
computed for a fixed chemical potential,this explains the jump of the current at / = ±0.4. 
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FIG. 3: The shift of the persistent current for the Majorana energy(two rings) cq = 0.01 . 

the result obtained here it is suggested that for two different fluxes, /i ^ the effect of 
the fluxes will be local . The correlation between the currents the two rings is negligible. 

In the past only measurements of few rings was possible . In recent years some of the 
experimental groups have claimed to measure the persistent current in a single ring. Some 
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FIG. 4: The shift of the persistent current for the Majorana energy (two rings) eo = 0.1 . 
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of the complications are related to the fact that the measurements have been done in diffusive 
limit (with 10-20 rings) and not in the ballistic limit which is considered in our calculation. 
The effect of the diffusive limit can be study by coupling the rings to a noise bath. 

XIII- Conclusions 

In conclusion the method of curved spaces introduced in the context of the Spin Hall effect 
has been applied to Topological Insulator and Weyl Semimetals. The method consists on the 
fact that in momentum space the coordinates are given by a momentum derivative, due to the 
fact that the spinor changes in the B.Z. the coordinate becomes covariant. The appearance 
of the spin connections in the covariant derivative generates the curvature. The Topological 
invariants, the Chern number are obtained after we impose the constraints time reversal, or 
charge conjugation ( for Superconductors). The Chern number are given by the covariant 
coordinates commutators. The formalism allows to obtain new Heisenberg equation of mo¬ 
tion which depends on the momentum curvature and real space curvature (magnetic fields). 
In particular the formalism is well suited to discuss T.I. or Weyl Semimetals on curved 
surfaces. The Weyl Semimetals are characterized by monopoles and antimonopoles which 
can be investigated using a new set of complete eigenfunctions coined harmonic monopoles. 
For the Topological Superconductors we considered the effect of the Majorana fermions on 
coupled metallic rings and computed the differential conductivity for the Andreev crossed 
reflection. 

Appendix A 

The physics of electrons in a periodic crystal is determined by the eigenvectors (spinors) 
\Un{k)) [n is the band-spin index) behavior in the Brillouin Zone k E (torus in a d 
dimensional momentum space). This behavior is similar to the parallel transport of a vector 
around a curve. We need to hnd the way the eigenvectors change under transport in the 
Brillouin Zone [121 EOl 1861 - 188] . The topological properties are encoded into the connection 
Ai (the vector potential in the momentum space) which measures the changes of \Un{k)) 
when it is transported in the Brillouin Zone. The changes are given by: 
id\Un{k)) — V'^^dk^\Um{k^)) = 0 (an index which appears twice implies a summation ). The 
matrix T™^ is given by = i{Un{k)\dki\Um{k)) = Af''^\k) where A is the connection. 
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Applying twice the (exterior) derivative we define the curvature = d{d\Un{k)) ( see 
Eqs.7.145 — 7.145 , Nakahara (2008) page 285 [ 13 ]) and find : 


F = d{d\Un{k))) 


+ r” (hrU(hl*' A dki \u,{k)) 


l[Fi,j]n,idk^ f\dk^\Ui(k)) 


the symbol A represents the wedge product 

{n,m)/J^\ -r 4(72,/)/ 


[F.,Um = d,A^;’^\k) - d,At^’^\k) + t[At^’\k),Af"^\k)], 
where Fij is the matrix curvature with the matrix elements [Eijjn.m 
commutator of the covariant derivative Ri = fi + , r, = id^i 


given in terms of the 

; [Ri, Rj] = Fij. 


Appendix B 


Topological invariant in two space dimensions nsing the topological invariant in 
four space dimensions. 


The topological response for time reversal invariant systems in one and two space di¬ 
mensions is not entirely clear . In three space dimensions we can use the Chern-Simons 
form A 3 (A, F)] to relate the the second Chern number C 2 in four space dimensions to three 
dimensions using the relation ch 2 = dlK^lA, F)]. In four dimensional momentum space 
the second Chern number C 2 is given by an index operator. In analogy with the index 
operator for the Dirac equation I introduce the index operator in the momentum space 
Ind. 


Ind. 


= Tr 


. -efd‘^kC{k)[iR4(k)\ C(k) 

ye 


|e -!>0 


Co 


(73) 


The operator zM is defined in terms of the non-Abelian spin connection : 


A‘/{k) = mk)\id^\Uf{k)) 

!K4(j:) = i ^ + A^(k)) = i'iaX‘‘{k) 

a=l,2,3,4 

5 T-i 0 

7® = 

0 1 

(74) 


7 ^ separates the conduction band from the valence band. 

In order to show that the index operator in four space dimensions is related to the index 
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operator in d=2 space dimensions, we introduce the transformation : 

^3(6) = g-\e,k)(^Aik) + d^g{e,k) (75) 

where 6 parametrizes the circle (see hgure 13.4 page 521 m)- Next we construct a family 
of gauge helds: 

= 0<f<l (76) 

The parameters {t, 9) form a disc with dD^ = . We construct fromU^ x 8“^ a manifold 

X S‘^ . We will call the patch [t, 9) the northern hemisphere Un and (s, 9) the south 
hemisphere Ug and the equator of corresponds to f = s = 1 (see the hgure with the 
two half sphere, hgure 13.4 page 521 in |I3]). The gauge potential in 2 + 2 dimensions can 
be written as: 


A]sf{t,9,k)= 0,0,A^'^{k)+g ^{9,k)deg{9,k) 

Ap^(t,9,k) = 0,0, ^(fc) 


(77) 


A{k) is the non-Abelian spin connection in two space dimensions. On the equator t = s = 1 
we have : 

— 9 ^ + d + (ig + dt)g 

where d is the exterior derivative [13] in two dimensions, de is the external derivative on 
and dtg = 0. A = (AAr,A 5 ) is the spin connection in 2 + 2 dimensions. We use the 
relation. 


det 


fZ)(^5(®))(fc)l = [det(*(^(fc)) 


1 

2 ^iw{A,d) 


(79) 


|def[zZl]| is gauge invariant and only det[iD] might have anomalous behavior . On the 
boundary disc dD^ = the phase dehnes the maping dD^ —)• 


tD{A)=t Yl + tA4k)P+), P+ = hl+7') 

a=l,2,3,4 


(80) 


\det[iD{A)] \ is gauge invariant and only the phase w{A,9) is anomalous and gives the wind¬ 
ing number ui, on the disc there are points at which det[iD{Ay\ vanishes. The index 


Ind. 


^ 2+2 


is given by. 


Ind. 


^ 2+2 


'S 2 xS 2 


C2(F) = Z/i 


(81) 


32 














where the curvature F is given by, 


F — (ci + (ig + dt)A. + 


(82) 


ui is the winding number which is is even or odd and corresponds to the index Z 2 introduced 
earlier. 

Following the procedure used before which relates the Chern character to the Chern-Simons 
form C 2 = d[K 3 {A, F)] we hnd: 


/ C2(F) — / C2(FAr)+ / C2(F5) — / 

/s2xS2 Jd ^ xs ^ Jd ^ xs ^ 751x52 


Since -^ 2 +i(As,F 5 ) = 0 we hnd: 

Ind. 


K2+i{^N, —-A2+i(A5, F5)|5=i 

(83) 


^ 2+2 


f iF 2 +i(A^,F^)= / K2+M^^^'^+9~^dgg,F3^^'^) = -{ — YTr[WdA] 

51x52 751x52 2 ! 7 r 

(84) 


where W = g ^{9)d0g{9). For 6* = 0 we hnd hF = 1 and Ind i]R 2+2 = is the winding 

number ) which is identical to the Z 2 index introduced by m- 

The procedure presented here allows a direct construction of the Z 2 invariant in two dimen¬ 
sions as an emergent object from four dimensions and therefore is related to the electromag¬ 
netic or sound wave response dehned in four space dimensions. Contrary to early procedures 
which used dimensional reduction the procedure proposed is based on deforming the spin 
connection to a family of higher dimensions gauge potentials 
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